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Abstract 

This article presents a technical innovation, study of solar power system based on the Stirling dish (SD) technology and design 
considerations to be taken in designing of a mean temperature differential Stirling engine for solar application. The target power source 
will be solar dish/Stirling with average concentration ratio, which will supply a constant source temperature of 320 °C. Hence, the system 
design is based on a temperature difference of 300 °C, assuming that the sink is kept at 20 °C. During the preliminary design stage, the 
critical parameters of the engine design are determined according to the dynamic model with losses energy and pressure drop in heat 
exchangers was used during the design optimisation stage in order to establish a complete analytical model for the engine. The heat 
exchangers are designed to be of high effectiveness and low pressure-drop. Upon optimisation, for given value of difference temperature, 
operating frequency and dead volume there is a definite optimal value of swept volume at which the power is a maximum. The optimal 
swept volume of 75 cm for operating frequency 75 Hz with the power is 250 W and the dead volume is of 370 cm . 

© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 

The harmony between environmental protection and 
economic growth has become a worldwide concern; there is 
an urgent need to effectively reuse solar energy, this source 
of energy is one of the more attractive renewable energy 
that can be used as an input energy source for heat engines. 
In fact, any heat energy source can be used with the Stirling 
engine. The solar radiation can be focused onto the heater 
of Stirling engine as shown in Fig. 1(a), thereby creating a 
solar-powered prime mover. The direct conversion of solar 
power into mechanical power reduces both the cost and 
complexity of the prime mover. In theory, the principal 
advantages of Stirling engines are their use of an external 
heat source and their high efficiency. Stirling engines are 
able to use solar energy that is a cheap source of energy. 

Studies about high temperature Stirling engines have been 
extensively reported in the literature [1] and commercial 
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units have been in operation for many years. On the other 
hand, low temperature Stirling engines are not as successful 
as their high temperature counterparts. However, the 
former have gained popularity in the last few decades due 
to this potential to tap a variety of low concentration 
energy sources available, such as solar. The increasing 
interest in Stirling engines is largely due to the fact the 
engine is more environmentally friendly than the widely 
used internal combustion engine, and also to its non¬ 
explosive nature in converting energy into mechanical form 
and thus leading to silent and cleaner operation, which are 
essential for special applications, such as military opera¬ 
tions and medical uses. 

The systems with very strong concentration [2] call upon 
an advanced and heavy technology, therefore are very 
expensive as they present, on the energy point of view, a 
limited interest. On the other hand, the systems without 
concentration are not economically viable. The best systems 
is with average concentration, leading to levels of temperature 
about 25CM50 °C, but very few work seem to be devoted to 
the installations with average concentration. The company 
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Nomenclature 

Subscripts 

A 

area, m 2 

c 

compression space 

Cp 

specific heat at constant pressure, Jkg -1 K -1 

ch 

load 

C pr 

heat capacity of each cell matrix, J K -1 

d 

expansion space 

C v 

specific heat at constant volume, Jkg -1 K -1 

E 

entered 

d 

hydraulic diameter, m 

ext 

outside 

D 

diameter, m 

f 

cooler 

dm 

wire diameter, m 

h 

heater 

fr 

friction factor 

moy 

mean 

Freq 

operating frequency, Hz 

P 

loss 

h 

convection heat transfer coefficient, 

Pa 

wall 


J m -2 s -1 K _1 

pis 

piston 

J 

annular gap between displacer and cylinder, m 

r 

regenerator 

k 

thermal conductivity, W m -1 K -1 

rl 

regenerator cell 1 

L 

length, m 

r2 

regenerator cell 2 

M 

mass of working gas in the engine, kg 

s 

left 

m 

mass flow rate, kgs -1 



m 

mass of gas in different component, kg 

Greek symbols 

NTU 

number of heat transfer unit 



P 

pressure, Pa 

e 

crank angle, rad 

Q 

heat, J 

8 

effectiveness 

Q 

power, W 

p 

Working GAS dynamic viscosity, kgm -1 s -1 

R 

gas constant, JkgK -1 

P 

density, kgm -3 

T 

temperature, K 

CD 

angular frequency, rads -1 

V 

volume, m 


mesh porosity 

W 

work, J 




BSR Solar Technologies GmbH, which developed the 
SUNPULSE, also works on a system intended to produce 
electricity starting from solar energy fairly concentrated, 
which leads to levels of temperature about 450 °C. 

Several analyses and simulation methods of the engine 
have been established [3], as well as the procedures for 
optimal design [4]. Most of the engines are fuel-fired and 
operate at high temperature, which highlights the need for 
careful material selection as well as good cooling system. 
For silent, light and portable equipment for leisure and 
domestic uses, low power engines may be more appropriate. 
Nevertheless, research in Stirling engine technology has been 
heavily masked by extensive and successful development of 
internal combustion engines, which have made Stirling 
engines less competitive. Hence, in order to design a low 
power engine using solar, new design specifications and 
optimisation criteria must be established [5-9]. This paper 
presents design considerations which may be taken to 
develop a solar Stirling engine with average concentration 
operating on mean temperature difference of 300 °C. 

2. Losses in a Stirling engine 

The energy losses in a Stirling engine are due to the 
thermodynamic and the mechanical processes. Compres¬ 
sion and expansion are not adiabatic. The exchangers are 
not ideal since the pressure drops in the engine and the 


losses of heat in the exchangers exist. To accurately predict 
power and efficiency requires an understanding of the 
principle parasitic loss mechanisms. 


2.1. Energy dissipation by pressure drops in heat exchangers 
dQvch 


Pressure drops due to friction and to area changes in 
heat exchangers is given by [10] 



2f^GV 
Ad 2 p 



0 1 * 

where G is working gas mass flow (kgm - s - ), d is the 
hydraulic diameter, p is gas density (kgm - ), V is volume 
(m ) and / r is the Reynolds friction factor. 

The internal heat generation which occurs when the gas 
is forced to flow against the frictional drag force, is given 
by [10]: 




m is the mass flow rate (kgs -1 ). 

The total heat generated by pressure drop in the different 
exchangers is 


^^Pch T d^Pchf T d^Pchrl T ^^Pchr2 "F ^^Pchh' 


( 3 ) 
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2.2. Energy lost by the internal conduction d()p cd 

Energy lost due to the internal thermal conductivity 
between the hot parts and the cold parts of the engine 
through the exchangers are taken into account. These 
losses are directly proportional to the temperature differ¬ 
ence at the ends of the exchanger; they are given for the 
different exchangers [11]: 


d Q 


Pcdh 



T r - h), 



k cd (Wm _1 K -1 ) is the material thermal conductivity; A is 
the effective area for conduction. 

So the total conduction loss is: 


dQpcdT — d£?Pcdr + ^0Pcdf + ^2pcdh* 



^2pcdr — ^cdr 





A f 

u 


(Tr_h — T f_ r ), 

(4) 

(Tf-r - r c _f), 

(5) 


2.3. Energy lost by external conduction d()p ext 

Energy lost by external conduction is considered in the 
regenerator which is not adiabatic. These losses are 
specified by the regenerator adiabatic coefficient, s ^ 1, 
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definite as the report between the heat given up in the 
regenerator by the working gas at its passage towards the 
compression space and the heat received in the regenerator 
by the working gas at its passage towards the expansion 
space [10]. So the energy stored by the regenerator at the 
time of the passage of gas from the expansion space to the 
compression space is not completely restored with this gas 
at the time of its return. 

For the ideal case of the regenerator perfected insulation, 
e = 1. 

The energy lost by external conduction is 

dgpext = (1 - <0(d&i + d Q t2 ). (8) 

The effectiveness of the regenerator s is given starting from 
the equation below [8] 

NTU 

£ “ 1 +NTU’ ( ^ 

NTU is the number of heat transfer unit: 


ntu = ^Ti, (10) 

Cpra 

where h is the overall heat transfer coefficient (hot stream/ 
matrix/cold stream), A wg refers to the wall/gas, or “wetted” 
area of the heat exchanger surface, C p the specific heat 
capacity at constant pressure, and m (kgs -1 ) the mass flow 
rate through the regenerator. 


2.4. Energy lost by Shuttle effect dg Pshtl 


Shuttling the displacer between hot and cold spaces 
within a machine introduces another mechanism for 
transferring heat from a hot to a cold space. Thus an 
important thermal effect appears in Stirling engines called 
‘Shuttle heat transfer’ having the effect of increasing the 
apparent thermal conductance loss. The displacer absorbs 
a quantity of heat from the hot source and restores it to the 
cold source. This loss of energy is given by [11]: 


dQpshtl — 


0.4Z fcpisAi (rjn r r ^ 

77 (T d q). 

JE d 



where J is the annular gap between displacer and cylinder 
(m), /c pis is the piston thermal conductivity (Wm -1 K -1 ), 
D d is the displacer diameter (m), L d is the displacer length 
(m), Z is the displacer stroke (m), T d and T c are, 
respectively, the temperature in the expansion space and 
in the compression space (K). 


3. Mathematical background 

There are many different ways to degrade the power 
produced by an ideal machine and to accurately predict 
power and efficiency requires an understanding of the 
design compartments. 

Mathematical model takes into consideration different 
losses and pressure drop in heat exchangers. 


Heat transfer and flow friction in the heat exchangers, 
i.e. the heater, the cooler and the regenerator, are evaluated 
using empirical equations under steady flow condition. 

No leakage is allowed either through the appendix gap 
or through the seals of the connecting rods. 

The temperature distribution in the various engine 
compartments is illustrated in Fig. 1(b). 

The gas temperature in the various engine compartments 
is variable. 

The cooler and the heater walls are maintained 
iso thermally at temperatures T paf and T pah . 

The pressure distribution is shown in Fig. 2. 

The gas temperature in the different compartments is 
calculated according to the perfect gas law: 



P c Vc 

Rm c 




PfVf 

Rm f 






PdVd 

Rm d 



The regenerator is divided into two cells rl and r2, each 
cell is been associated with its respective mixed mean gas 
temperature T r] and T r2 expressed as follows: 



Py 1 V rl 

Rm x \ ’ 



T ti — ~w~~ • (17) 

Rm r2 

An extrapolated linear curve is drawn through tempera¬ 
ture values T r \ and T r2 defining the regenerator interface 
temperature T r _ f , T r _ r and T r _ h , as follows [12]: 



3T r i — T y 2 



T r _ 


r—r 


T rl + Ty 2 




3T r 2 — T r i 
2 



According to the flow direction of the fluid, the interface’s 
temperatures: r c _ f , T f _ r , T r h and T h _ d are defined as 
follows [13]: 
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where r c _ f is the temperature of the interface between the 
compression space and the cooler, T f _ r is the temperature 
of the interface between the cooler and the regenerator, 
T r _ h is the temperature of the interface between the 
regenerator and the heater, T h _ d is the temperature of 
the interface between the heater and the expansion space. 
The matrix temperatures are so given by 

d^pari d O r | 


d t 


C pr d t 


d T 


par2 


dQ 


r2 


d t 


C pr d t 


( 21 ) 


( 22 ) 


where C pr is the heat capacity of each cell matrix (JK -1 ), 
Q r i is the quantity of heat exchanged to the regenerator rl 
(j), Qr 2 is the quantity of heat exchanged to the regenerator 
r2 (j), 7p ar i is the matrix temperature in the regenerator rl 
(K) and r Par2 is the matrix temperature in the regenerator 
r2 (K). 

By taking into account the conduction loss in the 
exchangers and the regenerator effectiveness, the power 
exchanged in the different exchangers is written 


dQf = ^f-4paf(Tpaf — Tf) — dg Pcdf . 


(23) 


d(2 r 2 — F7z r2 y 4par2(7~par2 F r2 ) 


d C?Pcdr2 


(24) 


dQ r i = Eh r \A par \(T VM \ — T t \) — 


d£?Pcdrl 


(25) 


nig, T e 


rh s , T§ 




■> 


Q 



m (0 P(t) 
T (t) v (t) 



w 


Fig. 3. Generalised cell. 


The work given by the cycle is 

dW „ dV c n dV d 

— — Pc -h P d — • 

dt dt dt 

The thermal efficiency given by the cycle is: 
W 




The total engine volume is: Vj = V c + Vf + V r \ + F r2 + 
V h + V d . 

The other variables of the dynamic model are given by 
energy and mass conservation equation, applied to a 
generalised cell as follows (Fig. 3): 


• Energy conservation equation : 


dQ h = /*h^pah(F p ah ~ Th) ~ d<2 Pcdh , (26) 

where d^pcdh is th e conduction loss in the cooler (W), 
d<2pcdri is the conduction loss in the regenerator rl (W), 
d£?pcdr 2 i s conduction loss in the regenerator r2 (W) and 
d^pcdh i s the conduction loss in the heater (W). 

The heat transfer coefficient of exchanges /zp, h r i, h r2 and 
h h is only available empirically [14]. 

The total exchanged heat is 

dQ = dQf + dQ vl + d() r2 + dQ h — dQ Pshtl . (27) 


de+C p r E m E -C p rsms = P— +C V ^df • ^ 3 °) 

Since there is a variable pressure distribution throughout 
the engine, we have arbitrarily chosen the compression 
space pressure P c as the baseline pressure. Thus, at each 
increment of the solution, P c will be evaluated from the 
relevant differential equation and the pressure distribution 
is determined with respect to P c . Thus it can be obtained 
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from the following expression: 


„ APf 

P{ ~ Pc+ 2 ’ 

(31) 

Ptl=P , + ^y^, 

(32) 

D D , (APrl + AP r2 ) 

Pr2 = P rl + 2 ’ 

(33) 

(AP rI + A A) 

Ph — P r2 T 

(34) 

p — p i 

Pd = Ph T 2 ‘ 

(35) 


Applying energy conservation equation to the different 
engine cells, we obtain: 

- C p 2Vf m cS = 1 (< C p P c ^ + C v V c , (36) 


dQ { - d Q pch{ + CpTc-f^fE - CpTt-rtkfs = 


C v V f d P c 


R d t 


(37) 


dO r | d(9p ch| .| CpCf-rWriE Cp T T — T fh T \ S — 


Cv V T \ d P c 


R At 


(38) 
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WrlS = 


C T 

^p J r—r 


dQrl ^^Pchrl ^p^f—r^rlE 


C v V r 1 dP c 


R d t 


(46) 


^r2S 
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C p T r _h 


(^dQ y2 ^Q?c\\v2 T 6^pT r _ r /h r 2E 


CyV r 2 dP c \ 
7? d? y’ 

(47) 


^hS 


1 

Cp TVd 


^dg h d2 Pchh + CpTr—hlWhE 


dmh 

dZfi 7 


C v F h dP c \ 
R dt y 

(48) 


where: m cS = mfe; m fS = m rm ; m r i S = m r2E ; m r2S = ^hE 
and mhs = ^dE- 


3.7. Solution method 

The systems of differential equations are written as 
follows: 

d Y = F(t,y), 

r(to) = ro, 

Y is a vector representing the unknown of each system, 
Y(t 0 ) = F 0 is the initial condition. 

These systems of equations are solved by the classical 
fourth-order Runge-Kutta method, cycle after cycle until 
steady. 



^Spchr2 T CpT r _ r m r 2E CpT r _E^2r2S — 


c v F r2 dP L 


p 


R 


dt 9 
(39) 


4. Design specification and concept 

4.1. Engine specification 


d Qh ^^Pchh T C p r r _ h m hE CpTh-e^hs — 


CyV h dP c 


R dr 


(40) 


C p 7V d m d - d<2 Pshtl =^(c p C 1 d C + c v F d ^ 


(41) 


Summing Eqs. (36)—(41) we obtain the pressure variation: 


d P c 


1 


dr Cy Ft 


7^(d Q — dQ FchT ) — 6 


p 


bW 

dr 


(42) 


• Mass conservation equation: 


M = md + m c + up + m r + n%. (43) 

The mass flow in the different engine compartments is 
given by the energy conservation Eqs. (36)—(41): 


^cs = - 


1 

RT c -f 




The engine parameters should be optimised [15] to avoid 
losses and to obtain high thermal efficiency for all the 
engine components especially heat exchangers. While the 
main target of the engine is to produce sufficient power to 
run a connecting application, there are conditions which 
pose critical constraints on the design, the working fluid is 
hydrogen and the temperature difference between the 
heater and the cooler is about 300 °C only. 

The engine presented in Fig. 4 uses a conventional crank 
mechanism driving two pistons by means of yoke linkage. 
The major feature of this is that there is almost no lateral 
movement of the connecting rods resulting in very small 
side forces on the pistons. With the lack of lateral 
movement of the connecting rods, there are relatively large 
unbalanced lateral forces due to the crankshaft counter¬ 
weight. Ross has a patented gear mechanism which 
balances the lateral forces by splitting and counter-rotating 
the counterweight 

4.2. Design concept 


mfs 



(aQ { - dO Pchr + C P 7V,/«,-£ 


C v Vf AP C \ 

R dt )’ 
(45) 


The yoke drive mechanism does not produce sinusoidal 
volume variations and the exact piston displacement 
functions are extremely complex. The volume variations 








































I. Tlili et at / Renewable Energy 33 (2008) 1911-1921 


1917 


are derived from geometric considerations in Fig. 5 and 
Table 1. 

The main Design concepts are listed in Table 2. 


5. Design analyses 

5.1. Relationship for engine power, swept volumes and dead 
volumes 

The purpose of this simulation is to estimate the main 
volumes of the engine spaces in terms of swept volumes and 



Table 1 

Volumes variations 


Geometrical parameters b\ = sin (p — r cos 9 

be = \Jb\ — {r cos 9) 2 
X = r sin 9 + be 


Displacements 


Volume variations 


dV c 
d 9 


= A p r 


cos 9 + sin 9 



+ 


rsin dcos 9 
be 


Y c = r[sin 9 — cos 9(b2/b\)] + be 
Y e = r[sin 9 + cos 9(b 2 /b\)] + be 

v c = V mc + T P (T max — r c ) 

V e = V me + A d (Y max - Ye) 


dV e 
d 9 


A d r 


r, ■ n {b2\ rsin 9 cos 9 

cos 9 — sin 9 ( —— I _|--- 

\b\ b e 


Table 2 

Concepts and target performance 


Parameters 

Values/type 

Engine type 

Alpha 

Working fluid 

Hydrogen 

Crank length 

r = 7.6 mm 

Yoke crank length 

b\ = 29 mm 

Piston length 

b 2 = 29 mm 

Displacement extremities 

F min = 17.75 mm 
F max = 39.28 mm 

Mean phase angle advance 

o 

O 

o\ 

II 

8 

Mass of gas in engine 

M = 0.35 g 

Hot space temperature 

T h = 590 K 

Cold space temperature 

T k = 290K 

Frequency 

Freq = 41.72 Hz 



x 


Fig. 5. Geometric derivation of the Ross Yoke drive equation. 









































































































































































































1918 


I. Tlili et al. / Renewable Energy 33 (2008) 1911-1921 


dead volumes using Dynamic model with losses since 
these factors are essential in estimating the preliminary 
configuration of the engine and will influence the sub¬ 
sequent optimisation process. Since it has been decided to 
adopt the successive alpha-type Ross Yoke configuration, 
the compression swept volume V c should be equal to the 
expansion swept volume V, d, and thus the swept volume 
ratio k = VJV d . 

In addition, at this stage, it is assumed that the mean 
pressure of the engine during operation is of 8.7 bar, which 
is the kind of pressure which normally occurs before the 
engine start-up. It is obvious from dynamic model 
equations that the net cycle power and the thermal load 
on the heat exchangers are direct linear functions of the 
engine speed (Operating rotation), the maximum pressure 
of the working fluid and the size of the engine, which is 
expressed in term of the swept volume [16]. However, the 
direct effects of the dead volume and swept volume to the 
engine power should be detailed. Figs. 6 and 7 illustrates 
the variation of the power as a function of the swept 
volume, which was calculated on the dead volumes of 535 
and 370 cm 3 under the fixed temperature difference 
of 300 °C. It is shown that the power increases when the 
swept volume increases until an optimal value. Also, it is 
noticeable that the power increases with the increase in 
speed. These two remarks imply that we have an optimal 
value of swept volume for maximum engine power for 
several speeds. By comparing the two graphs in Figs. 6 and 
7, based on the same swept volume, it can be said that the 
decrease in dead volume will lead to an increase in engine 
power. 

To illustrate the effect of the dead volume clearly, the 
variation of the engine power as a function of dead volume 
is calculated and the results are as shown in Figs. 8 and 9 
for the operating frequency of 75 and 35 Hz. 



Fig. 6. Relationship between swept volume and engine power (dead 
volume 535 cm ). 



Fig. 7. Relationship between swept volume and engine power (dead 
volume 370 cm ). 



Dead volume (cm 3 ) 

Fig. 8. Relationship between dead volume and engine power (frequen¬ 
cy = 75 Hz). 

From Figs. 6-9, it can be seen that the increase in the 
dead volume produces an exponential drop in the net 
power, which in turn decreases the maximum pressure. 
However, the calculation is performed under the assump¬ 
tion that the temperature difference is 300 °C, which can be 
obtained from the solar system with average concentration. 

5.2. Relationship for heater and cooler parameters 

An important factor in heat exchanger design is volume. 
Cooler and heater volumes contribute to large portions of 
dead volume. Previous studies showed that the dead 
volumes, which includes those in the heat exchangers, is 
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Dead volume (cm 3 ) 


Fig. 9. Relationship between dead volume and engine power (frequency = 
35 Hz). 

an essential factor in the Stirling engine design, where it 
should be small as possible [17]. To demonstrate the 
relationships for the heaters, specific conditions of 75 cm 3 
swept volume and 0.45 m tube length are used. 

After carrying out thermodynamic simulation for the 
heater, the variation of its tube diameter can be derived as 
a function of friction losses for several values of engine 
speeds as being depicted in Fig. 10 for the heater volume of 
165 cm . Similarly, Fig. 11 shows the graphs for the heater 
volume of 80 cm . Both graphs indicate an inverse 
proportionality between tube diameter and friction loss in 
the heater. The explanation of this variation is that the tube 
with smaller diameter having the same length delivers the 
same mass flux, thus generates a shorter entrance length 
and a thicker viscous boundary layer, which then leads to a 
higher friction factor of the flow. For the cooler, by using 
the same values for swept volume and tube length, the 
equivalent graphs for the cooler volumes of 165 and 80 cm 
are shown in Figs. 12 and 13, which indicates a similar 
pattern to that of the heater. 

In designing heat exchangers, an important considera¬ 
tion for the heat exchangers is to have an ability to supply 
or reject the required amount of heat to or from the engine. 
In this aspect, one crucial factor is the heat transfer area, 
which will decide the amount of heat energy to be 
transported. Hence, in order to achieve a high effectiveness 
for the heater and the cooler, larger transfer areas, and thus 
larger volumes, are needed. 

5.3. Relationship for regenerator parameters 

The effect of pressure drop in the regenerator of a mean 
temperature differential Stirling engine to thermal effi¬ 
ciency is very important since it can decrease the overall 
efficiency of the engine [16,17]. To analyse this effect and its 



0 0.5 1 1.5 2 2.5 3 

Tube diameter (cm) 

Fig. 10. Relationship of heater tube diameters with the friction losses 
(swept volume = 75 cm , tube length = 0.45 m, cooler volume 165 cm ). 



Fig. 11. Relationship of heater tube diameters with the friction losses 
(swept volume = 75 cm , tube length = 0.45 m, cooler volume 80 cm ). 

implication to the efficiency of the engine, six types of 
matrices has been selected and is being subjected to various 
pressure drops and engine speeds. The configurations for 
these six matrices are given in Table 3 for a standard total 
wire length of 5 m. The porosity of each matrix is 
important since it will have a direct impact on the 
performance of the regenerator, and can be determined 
by its geometry, namely, wire diameter, density of the mesh 
and the void volume. Any changes in the porosity will also 
change the regenerator effectiveness and the pressure drop, 
which eventually affects the engine efficiency. Therefore, 
the best matrix for the regenerator should possess both 
high efficiency and low-pressure drop. 
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Fig. 12. Relationship of cooler tube diameters with the friction losses 
(swept volume = 75 cm , tube length = 0.45 m, heater volume 165 cm ). 



Fig. 13. Relationship of cooler tube diameters with the friction losses 
(swept volume = 75 cm , tube length = 0.45 m, heater volume 80 cm ). 


Fig. 14 shows the relationship between the operating 
frequency and the pressure drop for these matrices. The 
pressure drop is found to be proportional to the frequency 
since an increase in frequency increases the mass flux 
through the regenerator as well as the pressure magnitude 
up to the same proportion for the same matrices. On the 
other hand, the decrease in mesh porosity leads to the 
higher friction factor as well as increases the pressure drop. 
Hence, it can be said that M1 has a lowest pressure drop in 
comparison to the others at a same speed because its 
porosity is the highest. In order to obtain a higher porosity, 
and thus the lower pressure drop, the meshes should be 
made from small wire diameter and should be as coarse as 


Table 3 

Geometrical properties of wire mesh for regenerator 


Matrix 

Wire diameter (m) 

Porosity (\jj) 

Ml 

0.0035 

0.9122 

M2 

0.0050 

0.8359 

M3 

0.0065 

0.7508 

M4 

0.0070 

0.7221 

M5 

0.0080 

0.6655 

M6 

0.0090 

0.6112 



possible. However, the pressure drop in the regenerator 
alone is not sufficient in deciding the best regenerator 
without considering its heat transfer behaviour. But 
Table 4 shows the relationship between the thermal 
efficiency, power of the engine and matrix type. The best 
matrix should compromise between high effectiveness and 
low-pressure drop in order to obtain minimal losses in the 
regenerator, and in this case, M6 with the porosity of 
0.6112 and wire diameter 0.009 m has been chosen for the 
design. 

The decrease in mesh porosity leads to the higher friction 
factor as well as increases the pressure drop. Hence, it can 
be said that Ml has a lowest pressure drop in comparison 
to the others at a same frequency because its porosity is the 
highest. In order to obtain a higher porosity, and thus the 
lower pressure drop, the meshes should be made from small 
wire diameter and should be as coarse as possible. 
However, the pressure drop in the regenerator alone is 
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Table 4 

Effect of matrix on power and thermal efficiency 


Matrix 

Power (W) 

Thermal efficiency (%) 

Ml 

159.39 

10.79 

M2 

226.42 

22.41 

M3 

249.20 

33.90 

M4 

252.50 

37.32 

M5 

255.92 

43.32 

M6 

256.77 

48.11 


not sufficient in deciding the best regenerator without 
considering its heat transfer behaviour this is the case of 
M6 In spite of the higher pressure drop we have better 
power and thermal efficiency because we have better heat 
transfer. 

The regenerator effectiveness s can be manipulated by 
varying wire diameter and wire length, which in turn 
changes the “wetted” surface area. It can be represented in 
form of the relationship between the porosity or the 
number of transfer units (NTU) and the thermal heating 
efficiency of the engine. If the wetted surface area is 
large, the resulting porosity should be low, and this 
provides the air or the work fluid with a large contacting 
surface to achieve a high rate of heat transfer. Hence, the 
NTU, and thus s, are increased when the surface area 
increases. The effect of s on the thermal efficiency of the 
engine is that it represents the ability to reject the heat to 
the working gas when the gas exits through the heater and 
the ability to absorb the heat when the gas exits through 
the cooler. 

6. Conclusion 

In this paper, a number of technical considerations in 
designing a mean temperature differential Stirling engine 
have been proposed. These considerations have been 
established through the use of the dynamic model with 
losses energy and pressure drop in heat exchangers. As a 
result, the optimal configuration for the design can be 
summarised as follow. 

• This studies show clearly that, for given value of 
difference temperature, operating frequency and dead 
volume there is a definite optimal value of swept 
volume ratio at which the power is a maximum. In 
this paper, the optimal swept volume is 75 cm for 
frequency = 75 Hz. 

• Upon optimisation, operating frequency has to be 
limited between 35 and 75 Hz at a temperature 
difference of 300 °C, where the best value is 75 Hz. 

• For the regenerator, its porosity plays a significant role 
in controlling pressure drop of the regenerator, which 
can be manipulated by varying wire diameter and 


length. For this engine, the selected parameters are 
the wire diameter of 3.5 mm with a total length of 5m 
and a porosity of 0.9122 to have low pressure drop but 
in our case M6 give the best thermal efficiency of the 
engine. 

• The heat exchanger volumes should be evaluated by 
considering both the pressure drop and the thermal 
efficiency of the engine. In our case the optimal heat 
exchanger volume has been found to be 165 cm for both 
the cooler and the tube dimension is 0.011 m in diameter 
and 0.450 m in length. 
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